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We present results of a detailed quantum mechanical study of a gas 

of N noninteracting electrons confined to a circular boundary and sub- 
ject to homogeneous dc plus ac magnetic fields {B = B^c + Bacfit), with 
/(t + 27r/a;o) = f{t))- We earlier found a one-particle classical phase dia- 
gram of the (scaled) Larmor frequency ujc = omegac/wo vs e = Bac/Bdc that 
separates regular from chaotic regimes. We also showed that the quantum 
spectrum statistics changed from Poisson to Gaussian orthogonal ensembles 
in the transition from classically integrable to chaotic dynamics. Here we find 
that, as a function of and (e,a;c), there are clear quantum signatures in 
the magnetic response, when going from the single-particle classically regular 
to chaotic regimes. In the quasi-integrable regime the magnetization non- 
monotonically oscillates between diamagnetic and paramagnetic as a function 
of N . We quantitatively understand this behavior from a perturbation theory 
analysis. In the chaotic regime, however, we find that the magnetization os- 
cillates as a function of N but it is always diamagnetic. Equivalent results are 
also presented for the orbital currents. We also find that the time-averaged 
energy grows like N'^ in the quasi-integrable regime but changes to a linear N 
dependence in the chaotic regime. In contrast, the results with Bose statistics 
are akin to the single-particle case and thus different from the fermionic case. 
We also give an estimate of possible experimental parameters were our results 
may be seen in semiconductor quantum dot billiards. 
Pacs 05.45. +b, 03.65.-w, 72.20. Ht 
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I. INTRODUCTION 



There is a long history of studies of the magnetic response of an electron gas, confined to a 
finite boundary. Starting with Bohr and van Leeuwen [Q, to Landau's finite diamagnetism in 
the quantum regime . The problem is still of current theoretical an experimental interest 

in particular due to the realization that for most geometries of the confining boundary 
one can find classically chaotic behavior [^- ||TH[. Most previous studies of this problem have 
assumed that the external magnetic field is static and they have concentrated in calculating 
the static magnetic susceptibility, except for the dynamic magnetic field experimental work 
of Reulet et al. [|l| . 

In an earlier paper |20| (referred to as I hereafter), we investigated the classical dynam- 
ics and the quantum signatures of classical chaos, for one electron confined to a circular 
quantum dot structure. The dot was subjected to uniform d.c. {Bdc) plus a.c. {Bacf{t)), 
with periodic f(t) = f(t + 27r/u;o)) perpendicular magnetic fields. There, we established an 
approximate phase boundary in the parameter space spanned by (e = Bac/ Bdc,u}c = ^cl^o) 
that separates the classically regular from the chaotic regimes, where ujc is the Larmor fre- 
quency of the d.c. field. The phase diagram shown in Fig. |^, which we shall often use in our 
analysis here, separates the quasi-integrable from chaotic regimes. In I we established clear 
correspondences between the transitions in the classical behavior and their corresponding 
quantum signatures. From the statistical properties of the quasienergy spectrum of the 
one-period evolution operator, going from Poisson to Gaussian orthogonal ensemble, to the 
semiclassical phase space correspondences via the Husimi quasienergy eigenfunction distri- 
bution functions. 

In this paper we present a detailed quantum mechanical study of the zero temperature 
magnetic response of a noninteracting electron gas confined to a circular boundary and 
subject to the same combination of a d.c. plus a.c. magnetic fields. Here we are interested 
in considering the magnetic response of this model for an electron system that satisfies 
the Pauli exclusion principle. Another basic question, first addressed in this paper, is how 
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does the transition from regular to chaotic behavior in the classical case, where the particles 
are indistinguishable, affect their fermionic quantum nature. Most previous studies of the 
quantum manifestations of classical chaos have centered on one-particle problems. Here we 
only address the important particle-statistics many-particle problem, and leave for a future 
study the relevant effects of electronic interactions. As we show below there are indeed 
clear manifestations of the particle statistics, which are different if we are in the classically 
integrable regime from those where the system is chaotic. 

The organization of the rest of the paper is as follows: In section II we briefly reca- 
pitulate the main elements of the single particle model studied in I, together with expres- 
sions for the matrix elements of the operators needed in our analysis. Next we outline our 
method to calculate the matrix elements of multi-electron operators in a basis of properly 
(anti) symmetrized eigenfunctions. In section III we present our main results for the magneti- 
zation, orbital currents and energy. We calculated both the time evolution of such operators 
and their time-averages as a function A^" and the parameters (e, tJc). We also include a 
perturbative calculation, fully described in the Appendix, that quantitatively explains our 
numerical results for the magnetization in the quasi-integrable weak-fleld regime. Finally, in 
section IV we present a summary of our conclusions, with an estimate of a few experimental 
parameters that may give an idea of the regimes in frequency and fields where the transition 
between integrable and chaotic regimes discussed in this paper could be tested. 

II. THE MODEL 

A. One-Electron Wavefunction 

We start by recalling the main features of the single particle formalism, as explained 
in paper I, and next its extension to the N non- interacting electron problem. The model 
we consider here is that of electrons confined to a disk, and subject to a steady (Bdc) and 
a time-periodic {Bac) magnetic field. After scaling to appropriate dimensionless units, the 



3 



model Hamiltonian considered here is, 

H = Hdc + Hi{T), (1) 

which in polar coordinates reads, 



and with the time- dependent kick component 



2 

n=—oo 



The dimensionless units are defined as. 



cD = — h- ^ e = — 



Here Rq is the radius of the disk quantum dot assumed to have rigid walls. To is the drive 
period of the a.c. field, uJc = e*Bdc/{m*c) is the static Larmor frequency, in terms of the 
effective electron mass m* (~ 0.067me), the screened electronic charge e* (~ 0.3e) ||2l|, and 
the dynamic frequency Uac = e*Bac/ {jn*c). 

The exact eigenf unctions of the static Hamiltonian Hdc are given in terms of the Whit- 
taker M functions [131, 



with n the principal quantum number, d. the angular momentum eigenvalue, and Nni a 
normalization constant. The frustration parameter, /, that measures the number of fiux 
quanta in the disk, is defined by, 

^ $0 " {hc/2e*) ~ h ^is' ' ^ ^ 

with is = (^^)^^^ the magnetic length and $o = ^ the quantum of fiux. The eigenenergies 

Ke = '2{xni + i), (8) 



are determined by the requirement that the wavefunction vanishes at the boundary i.e., by 
the zeros of the Whittaker function, M^^^ 1^1/2(2) = 0. 

We calculated the energy eigenvalues Eni for the static problem in a basis of Whittaker 
functions as a function of B^c, and checked our numbers by fully reproducing the results of 
Ref. 1^. The Whittaker functions have the advantage of being valid over the entire range of 
parameters, however, they are numerically difficult to evaluate for the full time- dependent 
problem. For convenience when calculating the time-dependent problem, we decided also to 
expand the total (single particle) wavefunction in a Fourier Sine-basis. In this case 

{r\Hdc |^n£(0)) = Enitpn£{r)^=, (9) 

V 2tt 



oo oo 

e 



n=l £=—oo 

^„,(r = l) = 0, [ ^Pl{r)rdr = l, (11) 

Jo 

[2 

and (r|'?/'„£) = Y - sin(n7rr). (12) 

This basis set is properly orthonormalized, and automatically satisfies the boundary 
conditions. To calculate the spectrum of the static problem, we used, nonetheless, the exact 
eigenvalues of Hdc, given by the zeros of the Whittaker functions. Doing this allowed us also 
to check the reliability of our Sine-basis numerical method. 

We then computed the required matrix elements of the operators we are interested in, 
within the Sine-basis method. For example, for the magnetization operator 

where C is the angular momentum operator and A(r) is the electromagnetic vector potential 
in normalized coordinates. In the present case, we take the magnetic field perpendicular to 
the plane, then the z-component of the magnetization operator is 
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with Lz the z-component of the angular momentum, fj,B — the Bohr magneton, and 
Mz the magnetization operator along the z-axis. The matrix elements of in the Fourier 
Sine-basis are given by. 



(f 8mn ^ 

-\o 2 ^(\^-"^n>- (15) 

Similarly, starting from the definition of the current density operator 

(where c.c. stands for complex conjugate), we have the following expression for the azimuthal 
current densities: — J^^"'""^ + where 

.(para) _ J^\_^ . „ „ „„J _ ^„ /.yN 

are the paramagnetic and diamagnetic current densities, respectively. In the Fourier Sine- 
basis, the matrix elements of the current densities (in units of K) are given by 

/ i7(p«-)i\ B~u I -Ci[2n7r] + 7E + ln(2n7r) (m = n) 

-Ci [(m -|- n)7r] -|- Ci [(m — n)7r] {m^ n), 

t.r \ f ~ \ h (m — n) 

2 [ ; (J-.), (m ^ n) 

where 7^ = 0.57721 566649 ... is Euler's gamma number, and Ci(x) is the Cosine integral. 

Finally, the expression for the one-period time-evolution operator C/^(t, Tq), for the single 
particle Hamiltonian, that satisfies the dynamical equation 

th^Ui{T,To) = iHdc + Hi{T))Ue{T,To), (20) 



is 



^7,(1,0) = exp ( -I ^r]A exp (-^Hdc ) • (21) 



The total (single particle) wavefunction at any integer multiple Nt of the period (hereafter 
taken to be 1), is given by repeated applications of Ui to the initial wavefunction: 



vl/,(r,0,iVr)) =f/f-|vl/,(r,0,O)). (22) 



B. Many-Electron Wavefunctions 



One can directly generalize the above single-electron formalism to the many-electron 
case. Take the initial A^-electron wavefunction to be 



\^r^,r2...,r^)) = \^l,2...,N)), 



(23) 



which is antisymmetric under exchange of an odd number of particles (the Pauli exclusion 
principle) : 



mi,...z,...j,...,N)) = -mi,...j,...z,...,N)). 



Let the i-th single-particle eigenstate satisfy the equation 



(24) 



(25) 



We know that the Slater antisymmetrization procedure for the non-interacting iV-electron 
state can be written as the following tensor product ||22|| : 



A • |<l>(l,2...,iV)) 



|^„,,,(i)) ®|^„,,,(i)) 

|^„,,,(2)) ®\^n,e,m 



®l^n^£iv(2)) 



I^PWj(l)) ® |^PW,}(2)) ■ ■ ■ ® |^p|„^,^}(iV)) 



(26) 



where A is the antisymmetrization operator, and P is the permutation operator 
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P{l,2,...,iV} = {Pl,P2,...,PiV}, (27) 

fl (even P) 

(28) 

-1 (oddP) 

The summation runs over all possible permutations. Furthermore, the trace of any sum of 
A^-body operators O = Xlili '^^^ ^e written as 



TT{d} = mi,...,N)\dmi,...,N)) 

1 ^ 

= ^Y.^'^-^Mm^n.im- (29) 



Now consider the time evolution of such a system. For simplicity, we take as the initial 
state the lowest energy (ground) state of the unperturbed system (i.e., without the a.c. field) 
allowed by the Pauli Principle : 

|$i,...,jv(t = 0)) = A ■ {|vl>„,,,(l)) ® \^n,i,{2)) ■ ■ • ® |vl>„,,,(iV))} . (30) 

The one-period time-evolution operator for the N-electron system is given by the tensor 
product, 

U = Ui{l, 0) ® f/2(l, 0) ■ ■ ■ ® Un{1, 0). (31) 

Since the antisymmetrization and time-evolution operators commute, the state after 
periods is simply given by 

|*l,...,7v(iVT)) = C/^"-|$l,...,7v(t = 0)) 

= A ■ {U^-\^^,,,il)) ■ ■ • ® U^-\^^,,,iN))}, (32) 

where we've used the notation f/j(l, 0) = Ui, for i = 1, ... , N. We can now generalize Eq. 
(p9|) for the trace of an operator at any integer multiple Nt of the period as, 

Tr{d}{t = Nt) = ($ i,...,7v(t = iVr)|6|<l>i,...,^(t = Nt)) 
1 ^ 

= ^ J2(^n,dJ,t = O)|(f/j)^-0,f/fn^„...(j,t = 0)). (33) 



In particular, for the average quantum time-dependent magnetization per electron we have 



'-"'jT^ = + E(*»,<.a)l(ci)"^(-('-|)r^f^|4'«.4(i)). (34) 

where L = YliLi ^i- Similarly, we can write the corresponding expressions for the time- 
dependent orbital currents, and the total time- dependent Hamiltonian average, which we 
term the averaged energy. For example, the time-averaged magnetization, ((M^)), is defined 
by 

((M,))= hm J2(M.){n). (35) 

^ n=l 



III. RESULTS 

We now come to the discussion of the main results of this paper, that are concerned 
with the dynamic and time-averaged properties of different relevant operators. The most 
striking features are observed in the magnetization of the system, which we shall discuss 
first as a function of the number N of electrons and (e, tUc). We also give results for the 
orbital current as well as interesting results for the time-averaged energy as a function of 
in different (e, uJc) parameter regimes. 



A. Magnetization, orbital currents and energy 

We start with the time-dependent dynamics of the magnetization and its corresponding 
power spectra for a single electron. The power spectrum S'(z/) is the square of the Fourier 
transform of the expectation value of the magnetization operator {Mz){t). (For notational 
simplicity, we write {M.^)(t) by (M)(t).) We keep lJc fixed, and sweep through values of e, 
from small to large. As can be seen from the phase diagram in Fig. |I], for fixed a)c, as e 
increases the underlying classical dynamics changes from quasi-integrable to chaotic. 

In Fig. ^(a) we see that (M) is diamagnetic in the regular regime and oscillates periodi- 
cally with time, reflected in the very strong peak in its power spectrum shown in Fig. |^(b). 
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As we increase the values of e (i.e., as we approach the chaos border in the phase diagram), 
the intermediate dynamics gets more complex, as shown in Figs, ^(a) and|^(b). We see in 
Fig. 0(b), that there are two peaks in S^u), which are due to the quasi-beats seen in Fig. 
0(a). In the chaotic region, < M > {t) shows essentially irregular behavior. Fig. ^(a), while 
S{i') has a broad background, shown in Fig. ^(b). Note that the average value of < M > (t) 
increases in magnitude, becoming steadily more diamagnetic in the chaotic regime. 

The one-electron behavior changes significantly with the addition of more electrons. The 
pattern of change from regular to chaotic is similar as in the one-electron case as we sweep 
through the same e-values as above. The spectral function develops more resonances in the 
regular region, whereas in the chaotic regime it has a broad band spectrum. As we continue 
to increase the number of electrons there are more "beats" in the time-dependence of the 
magnetization and more peaks in the spectral function. 

It is then more convenient to consider the time-averaged properties of the magnetization, 
the orbital current, or the total energy, as a function of the number of electrons. It is in 
this type of function that we can see important qualitative differences that represent the 
changes from the classical regular to chaotic behavior in the quantum dynamics. In Fig. 
^(a) we see that for two-electrons the time-averaged magnetization in the regular regime is 
paramagnetic, whereas for three electrons becomes diamagnetic again. We note that the dia- 
to para-magnetic changes are non-monotonic as a function of A^. For example, for A^ = 4 
it switches back to paramagnetic, but remains diamagnetic for both N = 5 and A^ = 6. A 
similar situation occurs with the orbital current as shown in Fig. ^(b). We mention that 
the specific value of the frustration parameter (/ = ^) determines if the magnetization 
flips from dia- to para-magnetic as we keep adding electrons. Basically, this phenomenon 
occurs when / ~ 0(1), i.e. when we add one flux quantum to the dot. In all other cases, 
the magnetization remains diamagnetic and monotonically increasing in magnitude as the 
number of electrons increases. We provide a theoretical perturbation theory explanation of 
these dia- to para-magnetic transitions result in the next subsection. 

As we increase the value of e, we enter the chaotic regime. There we find that for all elec- 
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tron numbers the magnetization is always diamagnetic, at least up to the maximum number 
of electrons we considered (~ 25). In Fig. ^(c) we show the time-averaged magnetization in 
the chaotic regime, which also oscillates as a function of A^, but it is always negative and of 
larger magnitude than in the quasi-regular regime. A similar situation occurs for the orbital 
current (as shown in Fig. |^(d), although it has less sharp changes as a function of than 
does << M ».) 

In Fig. ^(a) we consider the time- averaged magnetization for a fixed value of A^ = 1, 
e = 0.1 and = 0.1 as a function of Uc- In this quasi-integrable regime we see that << M » 
is diamagnetic and decays quadratically as a function of uJc- The situation changes in the 
chaotic regime, shown in Fig. |^(c), where there is also decay with uj^ but now the behavior 
is not as smooth as in the quasi-integrable regime. We show in Fig. ^(b) the behavior of 
the time averaged energy as a function of the number A^ of electrons. Here we see a clear 
quadratic growth as a function of A^. The situation is remarkably different when the single- 
particle classical dynamics is chaotic. In this case, shown in Fig. ^(d), the time-averaged 
energy grows clearly linearly with A^. This implies that the classically chaotic solutions do 
have a significant quantum signature in the averaged energy, that changes the quadratic 
quasi-integrable regime behavior to a linear A^ dependence in the chaotic regime. We now 
present a simple heuristic argument as to why the change over between quadratic and linear 
A^ behavior is actually directly related to the Pauli exclusion principle. We note that in the 
zero magnetic field case, each of the A^ electrons in the circular dot of radius Rq occupies 
an exclusion principle space of order Rq/N^^"^, while the static free particle kinetic energy 
changes like N/Rq. We expect that this situation does not change much when we are in 
the quasi-integrable regime, for finite fields and low frequencies. In the classically chaotic 
regime, in the presence of stronger magnetic fields or higher frequencies, the magnetic field 
will tend, on the average, to localize more the electrons to Larmor orbits inside the dot 
and in the boundaries. When the field is larger, so that the Larmor radius and Rq are 
comparable, the Landau levels have to be taken into account. In this case the electrons 
will not necessarily feel the presence of the boundary and they will remain localized in their 
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"chaotic" Landau orbits due to the time-dependent kicks. In this hmit the contribution from 
the kinetic energy is much less relevant, and the Larmor orbit radius will be less dependent 
of and Rq. 



B. Perturbative evaluation of the magnetization in the quasiregular regime 

In this subsection we present a perturbative analysis that provides an explanation for 
the magnetization oscillations as a function of in the quasi-integrable regime. Let us first 
consider the time-independent part of the one-electron Hamiltonian H , and write it as, 

Hdc = Ho + V, (36) 

where. 



and, 



V 4 (!) %=. (38) 

We will consider the limiting case of very small B^c field, i.e., <^ 1. As was first shown 
by Dingle, one can write the eigenvalues to first order in Uc, by considering the zero field 
basis functions of the disk Bessel eigenfunctions: 

^o|^l?)^i5i?|^S), (39) 
where the normalized eigenvalues are 0, 

^S'^f + (40) 

Here is the nth zero of the Bessel function Je{x), and the unperturbed basis functions 
are given by (we consider only the radial part, the angular part is clear) 

(r|V^i?)= J . Jiian,r). (41) 
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The matrix elements of the perturbation are then given by. 

2 I 1 



1 /cD, 



1 + ^-'^ 



m = n 



2V2; I 

The perturbed non degenerate eigen-functions are obtained from standard perturbation 
analysis, 



m^n mi nl 



IVff) + Y, J-p^^J'tC) + O (a') , (43) 



where we've used the unperturbed energy levels — E^^^ = ^{a'^i — ot^i)- Using Eq.( 



4 V'-*™.^ 

and the definition of the magnetization operator, Eg . ([l^) , the leading first order matrix 
element contribution to is given by 



2/ 



1 



2{P - V. 



+ 0{^t). (45) 



3 

Once again, if we take as the initial state the lowest energy state allowed by the Pauli 
Principle, 

|v|/i,...,^(0)) = A • {|^„,,,(1)) ® |^„,^,(2)) ■ ■ ■ ® \^n^,AN))} , (46) 



we can generalize Eq.(^) to the iV-electron case. We find that the averaged magnetization 
per electron, to first order approximation, is 

Next, we perform a linear-response theory analysis of the full time-dependent problem, 
assuming that e ^ 1. We show in the Appendix that within the perturbative approximation 
for a single electron, the average magnetization at time Nt is given by, 

(M,)(A^)^-|£, + ^(rVn. 
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where 



)n,ni — 2 



1 + 



2(£i 



'2 



+ 0{ui). (49) 



Clearly, the last term in Eq.(^) is of 0{u)^), so it can be ignored within the current approx- 
imation, which means that to lowest order, the time-dependence plays no significant role in 
determining the average magnetization. 

To test the approximation for {Mz)/N, Eg. (P7|) , we show in Fig. |^, a comparison of 
results from the perturbative and numerically exact calculations. The perturbative results 
agree remarkably well with the numerical calculations. We can now understand why the 
averaged magnetization oscillates in sign in the regular regime, and it is because {Mz)/N 
depends most strongly on Xlili^i- Whenever this sum of angular momentum quantum 
numbers fiips sign, so does the magnetization. For example, for the parameters shown in 
Fig. 0, the values of F{N) = ^i for successive values of are 

F{N = 1, 2, . . . , 12) = 0, -1, 0, -2, 0, 0, -3, 0, -1, 0, -4, 0, (50) 

corresponding to 

^, = (0, -1, 1, -2, 2, 0, -3, 3, -1, 1, -4, 4), (51) 

i.e., the values of the twelve-electron (unperturbed) ground state. This is a selection rule 
associated with the symmetries present in the system. 

What is interesting is that in the chaotic region, there is no such fiipping of the sign. This 
difference may constitute an experimentally accessible signature of chaos in the quantum 
system. Clearly, such behavior is exclusively a consequence of the Pauli principle, for we 
would not observe any change in the response per electron without it, since ignoring it in 
the multi-electron case would lead to a trivial rescaling of the single electron results. 



IV. CONCLUSIONS 



To summarize, we have studied a model of a non-interacting A^-electron system, con- 
fined to a circular structure with rigid boundaries, and subjected to perpendicular constant 

14 



and time-periodic magnetic fields. We studied the magnetization and orbital currents as 
a function of time, as well as the time-averaged magnetization ((M^)), and energy as a 
function of electron number A^. We can make a strong connection between the dynamic 
response of {Mz){t) (or it's power spectrum), and the underlying classical dynamics - as the 
classical system makes a transition to chaos as we vary the applied magnetic fields, the dy- 
namics changes from being harmonic to essentially noisy. There are three central significant 
conclusions: first, the Pauli Principle affects the behavior of this non-interacting system sig- 
nificantly, e.g. in terms of oscillations of ((M^)) as a function of A^. This behavior is directly 
related to the Pauli principle that allows the electrons to optimally reduce their averaged 
((M)) at specific values of the total angular momentum. Second, while these oscillations in 
the quasi-integrable regime cause the system to flip back and forth between dia- and para- 
magnetic behavior, the system remains diamagnetic at all times in the chaotic regime. We 
also found a very interesting change in the time-averaged energy as a function of A^, going 
from quadratic in the quasi-integrable regime to linear in the chaotic one. We provided a 
simple heuristic explanation of this behavior related to Pauli's exclusion principle. 

In this paper we have not considered the effects of Coulomb interactions that can signif- 
icantly complicate the analyses. There are static studies that have considered the changes 
in the classical dynamics due to interactions. What has been found in some examples is 
that if the system of non-interacting particles was non-chaotic, as the interaction parameter 
increases the dynamics can become chaotic p3|. In the quantum regime the Random Ma- 



trix Theory statistics can exhibit a transition from Poisson to orthogonal ensemble as the 



interaction strength increases . What happens in the time-dependent case considered in 
this paper that deals with quasi-energy statistics is not known at present. 

Here we have considered a circular disk in the presence of a time dependent magnetic 
field. It is only when we have the ac component of the field added to the dc one that chaos 
appears. In contrast, if the field is static but the geometry is changed one can have chaotic 
classical solutions. The relevance of the Pauli principle as seen in the zero temperature 
magnetization has been studied, for example, by At present we do not know what 
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happens when the geometry is not circular and we have a gas of Pauh electrons in the 
presence of a dc+ac magnetic field. We expect to consider the two problems mentioned 
above in the future. 

To conclude, we briefly give some estimates in terms of physical units of the field strengths 
and frequencies required to observe the effects predicted by our model calculations. In a 



GaAs-AlGaAs semiconductor the radius Rq of a quantum dot device [|10|, ||18[ can be between 
0.1 and 10/im, a sheet electron density n ~ 10^^ cm~^, a mobility /i ~ 265 000 cm^/V-s, and 
a characteristic level spacing Ae ~ 0.05 meV or ~ 0.5 K. In the ballistic electronic motion 
regime the elastic mean free path ~ lOyum, with phase coherence length varying between 
15 and 50 /im. Typically the power injected is smaller than 1 nW, which is necessary to 
avoid electron heating. For a dot radius of Rq ~ l/im, the kick frequency Uq can be obtained 
from Eqs.(ll) as ujq = — 1 — | GHz. Then the required B^c and Bac magnetic fields have 
the values: Bdc = ^^^^^^Cjc — 20^ gauss, and Bac = eBdc ^ 20^ gauss. The a.c. Larmor 
frequency is oOac = — 20 ^ MHz. With these values, in the quasi-integrable regime, 
with parameters (e, cDc)*-^*^^^ = (0.1,0.1), we get uJq^^^ ~ 20 GHz and Bac^^ — 20 gauss. 
In the chaotic regime we take the parameters {e,Uc)^'^'^°'°^^ = (2.0,2.0). which leads to 
^{chaos) ^ 20 QHz and B^c"""^^ — 800 gauss. These results for the regular and chaotic 
regimes are within experimental reach. 
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APPENDIX A: 

In this appendix we provide the derivation of Eq.(^). For any operator A, the expec- 
tation value of the hnear response under the action of a constant Hamiltonian Hq and a 
time-dependent perturbation V(t) is given by 



1 



(A)(t) = (A)o + - / dt,{[A{U-t),V{U)])o, (Al) 







where ( ■ )o = Tr{po-}, and po is the density matrix associated with the unperturbed Hamil- 
tonian Hq. Thus, 



{[A{t, - t), V{h)])o = Tr [po[A{h - t), V{h)] 

- poV'(ti)e'^°^*^"*^/''ie-'^" . (A2) 
For a single particle pure state lipm), Po = |^ni)('^nj- Thus, from above, 

n 



n 



In our case. 



n 



h 2 

Vih) = vY,5ih-p). (A4) 

p=0 



Since the {iV'n)} sue real and A and V are Hermitian, 
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where 



A 



f 



f 



ni,n 



K,n, = ^(^) {r^)n,n,^S{ti-p), and 



2 V 2 



ni,n 



En, — En 

h ■ 



Finally, we have 



(...)o = -^ 



X/n )n,ni ( ^i^ni,n + q ('^ )n,ni 

EfJo sin {^n.,n(il - t)} - p) 



(A5) 



(A6) 
(A7) 

(A8) 



(A9) 



and, 



{A)o={iJn^M,\^n,) = -£i 



— I 2\ 

2 V Irii 



(AlO) 



Substituting Eqs. (|A9|) and (|A10| ) into Eq.(|Al[), we get Eg. p8|) . This completes our deriva- 
tion of the linear response result. 
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FIG. 1. Classical phase diagram separating the regular from chaotic regions, used in the analysis 
of this paper. See text for the definition of the dimensionles variables in the axes. 

FIG. 2. (a) Single electron time-dependent magnetization (M)(t) for e = 0.1, ujc = 0.1, and 
h = {).l in the quasi-regular regime (see Fig. ^). (b) Power Spectrum, S{v), of (Af) corresponding 
to the parameter values of Fig. |2|(a). The scale of v frequencies is arbitrary. In all figures the axis 
are given in terms of dimensionless units. 

FIG. 3. (a) Same as in Fig. ^(a) for e = 1.0, cJc = 1.0, h = 0.1, which is in the regime that is 
approaching the chaotic region in parameter space (see Fig. |l|). (b) Same as in Fig. |2|(b) for the 
parameters of Fig. |^(a). 

FIG. 4. (a) The magnetization as a function of time for the parameters e = 2.0, cJc = 2.0, 
h = 0.1, deep in the chaotic regime (see Fig. ^). (b) Here the classical chaos is clearly revealed in 
the noisy structure of < M > and the broad band spectrum of S{v). 

FIG. 5. (a) Time-averaged magnetization ((M)) (per electron) in the quasi-integrable regime for 
e = 0.1, tOc = 0.1, h = 0.1, as a function of electron number A^. We see that ((M)) oscillates non 
monotonically with A^. See text for a theoretical perturbative explanation of these oscillations, (b) 
Total time-averaged orbital magnetization (per electron) as a function of N for the same parameter 
values as in Fig. 1(a). Here we see, as expected, a direct correspondence with the magnetization 
of Fig. ^(a). (c) Same as in Fig. |5|(a) for parameters in the chaotic regime, e = 2.0, Uc = 2.0, 
/i = 0.1. In this case ((M)) oscillates but its behavior is purely diamagnetic. (d) Orbital current for 
the chaotic parameter values of Fig. ^(c), with a clear correspondence to the purely diamagnetic 
nature of the magnetic response of the system in the chaotic regime. 
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FIG. 6. (a) Time-averaged magnetization for one electron as a function of Cbc for e = 0.1, and 
^ = 0.1, with increment in diamagnetism as Cbc increases, (b) Time-averaged energy as a function 
of electron number, with a clear quadratic dependence, for the same parameters as in (a), (c) 
((M)) in the chaotic regime for one electron, with e = 10, ^ = 0.1, as a function of ujc- Here we 
also note an increase of diamagnetism as (be grows, (d) {{E)) for the same parameters as in (b) in 
the chaotic regime with a clear linear dependence on A^. See text for further discussion of these 
results. 

FIG. 7. Comparison between exact (□) and perturbative (+) calculations for the time-averaged 
magnetization as a function of electron number. The parameters here are, e = 0.01, ^ = 0.1, and 
u)c = 0.01. Note the almost exact agreement between the two calculations, with the magnitude of 
the diamagnetic response much smaller than that of the paramagnetic response. 
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